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Centrality Dependence of Charged Particle Production at Large 
Transverse Momentum in Pb-Pb Collisions at ^snn = 2.76 TeV 



The inclusive transverse momentum (pj) distributions of primary charged particles are measured in 

the pseudo-rapidity range |tj | < 0.8 as a function of event centrality in Pb-Pb collisions at -^/snn = 2.76 TeV 



with ALICE at the LHC. The data are presented in the pj range 0.15 < p-\ < 50 GeV/c for nine cen- 
trality intervals from 70-80% to 0-5%. The Pb-Pb spectra are presented in terms of the nuclear mod- 
ification factor Raa using a pp reference spectrum measured at the same collision energy. We observe 
that the suppression of high-px particles strongly depends on event centrality. In central collisions 
(0-5%) the yield is most suppressed with Raa ~ 0. 1 3 at pj = 6-7 GeV/c. Above pj = 7 GeV/c, there 
is a significant rise in the nuclear modification factor, which reaches Raa ~ 0.4 for pj > 30 GeV/c. 
In peripheral collisions (70-80%), the suppression is weaker with Raa ~ 0.7 almost independently 
of pj. The measured nuclear modification factors are compared to other measurements and model 
calculations. 
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1 Introduction 

High-energy collisions of heavy-ions enable the study of hot and dense strongly interacting matter JT|-0. 
At sufficiently high temperature, it is expected that partons (quarks and gluons) are the dominant degrees 
of freedom. During the very early stage of the collision, some of the incoming partons experience scat- 
terings with large momentum transfers. These partons lose energy when they traverse the subsequently 
formed hot and dense medium. One of the major goals of the heavy-ion physics programme at the LHC 
is to understand the underlying mechanisms for parton energy loss and use this as a tool to probe the 
properties of the medium. 

Parton energy loss in heavy-ion collisions was first observed at RHIC as the suppression of high-px 
particle production in Au-Au collisions compared to expectations from an independent superposition of 
nucleon-nucleon collisions [6-9]. At RHIC, the particle production in central (0-5%) Au-Au collisions 
at ^snn = 200 GeV is suppressed by a factor of 5 at pj = 5-6 GeV/c JUS], an d is consistent with being 
independent of pj over the measured range 5 < pj < 20 GeV/c iflQl . 

The increase of the charged particle density (dAf c h/dT7) at mid-rapidity from RHIC energies to actual 
LHC energies by a factor of around 2.2 [11] implies a similar increase in energy density. However, the 
observed suppression of high-/?x particle production also depends on the ratio of quarks to gluons due 
to their different color factors, and on the steepness of the pj spectra of the scattered partons. At the 
LHC the initial parton pj spectra are less steep than at RHIC and the ratio of gluons to quarks at a given 
pj is higher lPT2l . The measurement of high-/?T hadron production at the LHC helps to disentangle the 
effects which cause the suppression and provides a critical test of existing energy loss calculations 021. 
In particular, the large pj reach provides a means to study the dependence of the energy loss on the initial 
parton energy. 

We present a measurement of the pj distributions of charged particles in 0.15 < pj < 50 GeV/c with 
pseudo-rapidity |tj| <0.8, where T] = —In [tan (6/2)}, with 6 the polar angle between the charged particle 
direction and the beam axis. Results are presented for different centrality intervals in Pb-Pb collisions 
at -y/snn = 2.76 TeV. They are compared with measurements in pp collisions, by calculating the nuclear 
modification factor 



d 2 < A /dTjd PT 

where A^ A and represent the charged particle yield in nucleus-nucleus (AA) collisions and the cross 
section in pp collisions, respectively. The nuclear overlap function Taa is calculated from the Glauber 
model Ifl4l and averaged over each centtality interval. 

Early results from ALICE |[T5l showed that the production of charged particles in central (0-5%) Pb-Pb 
collisions at ^Ann = 2-76 TeV is suppressed by more than a factor of 6 at p t = 6-7 GeV/c compared 
to an independent superposition of nucleon-nucleon collisions, and that the suppression is stronger than 
that observed at RHIC. The present data extend the study of high-px particle suppression in Pb-Pb out 
to pj = 50 GeV/c and include a systematic study of the centrality dependence. Moreover, the systematic 
uncertainties related to the pp reference were significantly reduced with respect to the previous measure- 
ment by using the pj distribution measured in pp collisions at ^Js = 2.76 TeV lfl6l . 



2 Experiment and Data Analysis 

The ALICE detector is described in IfTTll . The Inner Tracking System (ITS) and the Time Projection 
Chamber (TPC) are used for vertex finding and tracking. The minimum-bias interaction trigger was 
derived from signals from the forward scintillators (VZERO), and the two innermost layers of the ITS 
(Silicon Pixel Detector - SPD). The collision centrality is determined using the VZERO. In addition, 



2 



The ALICE Collaboration 



Table 1: Average values of the number of participating nucleons (N part ) and the nuclear overlap function (Taa) lfl4l 
for the centrality intervals used in the analysis. 



V^CllLldllL y 


VVpart/ 


\ 7 AA/ \ lnu ) 


0-5% 


383 ± 3 


26.4 ± 1.1 


5-10% 


330 ±5 


20.6 ± 0.9 


10-20% 


261 ±4 


14.4 ± 0.6 


20-30% 


186 ±4 


8.7 ± 0.4 


30-40% 


129 ± 3 


5.0 ± 0.2 


40-50% 


85 ±3 


2.68 ±0.14 


50-60% 


53 ±2 


1.32 ±0.09 


60-70% 


30.0 ± 1.3 


0.59 ± 0.04 


70-80% 


15.8 ± 0.6 


0.24 ± 0.03 



the information from two neutron Zero Degree Calorimeters (ZDCs) positioned at ±114 m from the 
interaction point was used to remove contributions from beam-gas and electromagnetic interactions. The 
trigger and centrality selection are described in more detail in [ 1 1 ]. 

The following analysis is based on 1.6- 10 7 minimum-bias Pb-Pb events recorded by ALICE in 2010. 
For this study, the events are divided into nine centrality intervals from the 70-80% to the 0-5% most 
central Pb-Pb collisions, expressed in percentage of the total hadronic cross section. The event central- 
ity can be related to the number of participating nucleons _/V part and the nuclear overlap function Taa by 
using simulations based on the Glauber model lfl4l . The average values of N pa[t and Taa for each cen- 
trality interval, (N pan ) and {Taa}, along with their corresponding systematic uncertainties, are listed in 
Table [T] The errors include the experimental uncertainties on the inelastic nucleon-nucleon cross section 
a ijiei = 64 ± 5 mb at a/snn = 2.76 TeV |[T8l and on the parameters of the nuclear density profile used in 
the Glauber simulations (more details in iflTlD . 

The primary vertex position was determined from the tracks reconstructed in the ITS and the TPC by 
using an analytic % 2 minimization method, applied after approximating each of the tracks by a straight 
line in the vicinity of their common origin. The event is accepted if the coordinate of the reconstructed 
vertex measured along the beam direction (z-axis) is within ±10 cm of the nominal position. The event 
vertex reconstruction is fully efficient for the event centralities covered. 

Primary charged particles are defined as all prompt particles produced in the collision, including decay 
products, except those from weak decays of strange hadrons. A set of standard cuts based on the number 
of space points and the quality of the momentum fit in the TPC and ITS is applied to the reconstructed 
tracks. Track candidates in the TPC are required to have hits in at least 120 (out of a maximum of 159) 
pad-rows and % 2 per point of the momentum fit smaller than 4. Such tracks are projected to the ITS 
and used for further analysis if at least 2 matching hits (out of a maximum of 6) in the ITS, including at 
least one in the SPD, are found. In addition, the % 2 per point of the momentum fit in the ITS must be 
smaller than 36. In order to improve the purity of primary track reconstruction at high pj we developed 
a procedure where we compare tracking information from the combined ITS and TPC track reconstruc- 
tion algorithm to that derived only from the TPC and constrained by the interaction vertex point. We 
calculated the Ztpc-its between these tracks using the following formula 

^TPC-ITS =( V TPC — VTPC ITS) T ' (CtPC + CxPC ITs) 1 ' (VXPC ~ VXPC ITs) (2) 

where vxpc, vjpc-its an d Cxpc, Cxpc-its represent the measured track parameter vectors v = (x,y,z, 0, 0, 1//?t) 
and their covariance matrices, respectively. If the ^| PC -its * s larger than 36 the track candidate is re- 
jected. At pi = 0.15-50 GeV/c, this procedure removes about 2-7% (1-3%) of the reconstructed tracks 
in the most central (peripheral) collisions. This procedure in fact removes high-px fake tracks, which 
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Table 2: Contribution to the systematic uncertainties on the pj spectra (0.15-50 GeV/c) for the most centtal and 
peripheral Pb-Pb collisions. Also listed are the systematic uncertainties on the pp reference (0.15-50 GeV/c) 1 16 1. 





0-5% 


70-80% 


Centrality selection 


0.5% 


2.6% 


Event selection 


3.2% 


3.4% 


Track selection 


4.1-7.3% 


3.6-6.0% 


Tracking efficiency 


5% 


5% 


pj resolution correction 


<1.8% 


<3% 


Material budget 


0.9-1.2% 


0.5-1.7% 


Particle composition 


0.6-10% 


0.5-7.7% 


MC generator 


2.5% 


1.5% 


Secondary particle rejection 


<1% 


<1% 


Total for pj spectra 


8.2-13.5% 


8.3-11.9% 


Total for pp reference 

pp reference normalization 


6.3-18.8% 
1.9% 



originate from spurious matches of low pj particles in the TPC to hits in the ITS, and would result in an 
incorrect momentum assignment. 

Finally, tracks are rejected from the sample if their distance of closest approach to the reconstructed 
vertex in the longitudinal direction d z is larger than 2 cm or dxy ^ 0.018cm -)- 0.035cm • in the 
transverse direction with pj in GeV/c, which corresponds to 7 standard deviations of the resolution in 
d xy (see |fT9l for details). The upper limit on the d z (d z < 2 cm) was set to minimize the contribution 
of tracks coming from pileup and beam-gas background events. These cuts reject less than 0.5% of the 
reconstructed tracks independently of pj and collision centrality. 

The efficiency and purity of the primary charged particle selection are estimated using a Monte Carlo 
simulation with HIJING EOl events and a GEANT3 [21] model of the detector response. We used a 
HIJING tune which reproduces approximately the measured charged particle density in central colli- 
sions [11]. In the most central events, the overall primary charged particle reconstruction efficiency 
(tracking efficiency and acceptance) in |tj| < 0.8 is 36% at pj = 0.15 GeV/c and increases to 65% for 
pj > 0.6 GeV/c. In the most peripheral events the efficiency is larger than that for the central events 
by about 1-3%. The contribution from secondary particles was estimated using the d xy distributions of 
data and HIJING and is consistent with the measured strangeness to charged particle ratio from the re- 
construction of K^, A and A invariant mass peaks in Pb-Pb Il22l . The total contribution from secondary 
tracks at pj = 0.15 GeV/c is 13 (7)% for central (peripheral) events and decreases to about 0.6% above 
pi = 4 GeV/c for both central and peripheral events. From a systematic variation of the #| PC _ IT s cut 
and comparison of track properties in MC to data we conclude that the number of properly reconstructed 
tracks rejected as high-pj fake tracks is around 1-2% (0.5-1%) in the most central (peripheral) colli- 
sions. We also conclude that the contribution from the high-pj fake tracks to the pj spectra is negligible 
independently of the collision centrality and pj. 

The transverse momentum of charged particles is reconstructed from the track curvature measured in the 
magnetic field B = 0.5 T using the ITS and TPC detectors. The pj resolution is estimated from the track 
residuals to the momentum fit and verified by cosmic muon events, and the width of the invariant mass 
peaks of A, A and K ( s ' reconstructed from their decays to two charged particles. For the selected tracks 
the relative pj resolution {o{pj)/pj) amounts to 3.5% at pj = 0.15 GeV/c, has a minimum of 1% at 
pj = 1 GeV/c, and increases linearly to 10% at pj = 50 GeV/c. It is independent of the centrality of 
the selected events. From the study of the invariant mass distributions of A and K° as a function of pj 
we estimate that the relative uncertainty on the pj resolution is around 20%. From the mass difference 
between A and A and the ratio of positively to negatively charged tracks, assuming charge symmetry 
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Fig. 1: Charged particle pj distribution measured in Pb-Pb collisions in different centrality intervals. The spectra 
are scaled for better visibility. The dashed lines show the pp reference [16| spectra scaled by the nuclear overlap 
function determined for each centrality interval (Table[T} and by the Pb-Pb spectra scaling factors. The systematic 
and statistical uncertainties for Pb-Pb are added quadratically. The uncertainties on the pp reference are not shown. 



at high pj, the upper limit of the systematic uncertainty of the momentum scale is estimated to be 
|A(/?t) / Pt\ < 0.005 at pi = 50 GeV/c. This has an effect of around 1.5% on the yield of the measured 
spectra at the highest pj. To account for the finite pi resolution, correction factors for the reconstructed 
pj spectra at pj > 10 GeV/c are derived using a folding procedure. The corrections depend on collision 
centrality due to the change of the spectral shape and reach 4 (8)% at p t = 50 GeV/c in the most central 
(peripheral) collisions. 

The systematic uncertainties on the pj spectra are summarized in Table |2] They were estimated by a 
comparison of different centrality measures using hits in the SPD instead of VZERO and by varying 
the event and track quality selection criteria and the Monte Carlo assumptions. In particular, we studied 
a variation of the most abundant charged particle species (pions, kaons, protons) by ±30% to match 
the measured ratios and their uncertainties |[22l . The material budget was varied by ±7% ||23], and the 
secondary yield from strangeness decays in the Monte Carlo by ±30% to match the measured d xy distri- 
butions. Moreover, we used a different event generator, DPMJET E4l . to calculate MC correction maps. 
The systematic uncertainties on the pj spectra, related to the high-/?T fake track rejection procedure, 
were estimated by varying the track matching criteria in the range 25 < Xtpc-its < 49, and amount to 
1-4% (1-2%) in the most central (peripheral) collisions. The total systematic uncertainties on the cor- 
rected pj spectra depend on pj and event centrality and amount to 8.2-13.5% (8.3-11.9%) in the most 
central (peripheral) collisions. 
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Fig. 2: Nuclear modification factor Raa of charged particles measured in Pb-Pb collisions in nine centrality in- 
tervals. The boxes around data points denote /^-dependent systematic uncertainties. The systematic uncertainties 
on the normalization which are related to (7Xa) and the normalization of the pp data are added in quadrature and 
shown as boxes at 7?aa = 1 ■ 

A dedicated run of the LHC to collect pp reference data at ^fs = 2.76 TeV took place in March 201 1. 
Data taken in this run were used to measure the charged particle pj spectrum that forms the basis of the 
pp reference spectrum for Raa- Using these data the systematic uncertainties in Raa related to the pp 
reference could be significantly improved (Tabled compared to the previous publication [15], allowing 
for an exploration of high-px particle suppression in Pb-Pb out to 50 GeV/c. More details about the pp 
reference determination can be found in lfl6ll . 



3 Results 

The fully corrected pj spectra of inclusive charged particles measured in Pb-Pb collisions at y/^m = 2.76 TeV 
in nine different centrality intervals, and the scaled pp reference spectra are shown in Fig.Q] In the most 
peripheral Pb-Pb collisions the shape of the pj distribution is similar to that observed in pp collisions. 
A marked depletion of the spectra at high transverse momentum (pj > 5 GeV/c) develops gradually as 
centrality increases, indicating strong suppression of high-px particle production in central collisions. 

The nuclear modification factors for nine centrality intervals are shown in Fig. [2] In the most peripheral 
collisions (70-80%), the suppression is independent of pi, with Raa ~ 0.7. In central collisions, a pro- 
nounced minimum at about pj = 6-7 GeV/c is observed, while above p-\ = l GeV/c there is a significant 
rise of the nuclear modification factor. This rise becomes gradually less steep with increasing pj. In the 
most central collisions (0-5%), the yield is most suppressed, Raa ~ 0.13 at pj = 6-7 GeV/c, and Raa 
reaches « 0.4 for p T > 30 GeV/c. 
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Fig. 3: Nuclear modification factor /?aa of charged particles as a function of (N part ) (top panel) and diVch/drj 
(bottom panel) measured by ALICE in Pb-Pb collisions in different px-intervals, compared to PHENIX results 
in 5 < pj < 7 GeV/c [9|. The boxes around the data represent the /^-dependent uncertainties on the Pb-Pb pj 
spectra. The boxes at R^a = 1 represent the systematic uncertainties on the pp reference in different /^-intervals 
(pT-interval increases from left to right, the left-most is for PHENIX). The systematic uncertainties on the overall 
normalization for ALICE and PHENIX are not shown. 

The dependence of Raa on the collision centrality, expressed in terms of N part and the charged particle 
multiplicity density (dA^h/dr/), are shown in Fig. [3] for different intervals of pj. Also shown are results 
from PHENIX at RHIC in Au-Au collisions at ^/snn = 200 GeV [9]. The strongest centrality depen- 
dence is observed for particles with 5 < pi < 7 GeV/c. At higher pj, the centrality dependence weakens 
gradually. In comparison to results from RHIC, the LHC data in the same pi window show a suppression 
which is larger by a factor of about 1.2 at all (N part ) (Fig. |3l top panel). This implies that the shape of 
the A^p art dependence at RHIC and the LHC is very similar when the same pi is compared, indicating a 
strong relation between collision geometry and energy loss. The overall increase of suppression at the 
LHC as compared to RHIC may be expected from the larger density and longer lifetime of the fireball. 
The suppression reaches similar values when results from RHIC are compared to results from the LHC 
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Fig. 4: Nuclear modification factor Rpj± of charged particles measured by ALICE in the most central Pb-Pb 
collisions (0-5%) in comparison to results from CMS [25 1 and model calculations [26-30]. The boxes around 
the data denote /^-dependent systematic uncertainties. For CMS statistical and systematic uncertainties on /?aa 
are added in quadrature. The systematic uncertainties on the normalization which are related to (Taa) and the 
normalization of the pp data are added in quadrature and shown as boxes at Raa = 1 (the right-most is for CMS). 

in terms of (oWch/dTj), as shown in Fig. [3] (bottom panel). Larger values of suppression than at RHIC 
are observed in central collisions at the LHC, where the charged particle multiplicity exceeds that of the 
most central collisions at RHIC. It should be noted that the suppression at a given centrality results from 
a subtle interplay between the parton pj spectrum, the quark-to-gluon ratio, and the medium density, 
all of which exhibit a significant energy dependence. Further model studies are needed to evaluate their 
relative contributions. 

The ALICE measurement of Raa in the most central Pb-Pb collisions (0-5%) is compared to the CMS 
result [25 ] in Fig. [4] Both measurements agree within their respective statistical and systematic uncer- 
tainties. 

In Fig. [4j the measured Raa for 0-5% central collisions is also compared to model calculations. All se- 
lected models use RHIC data to calibrate the medium density and were available before the preliminary 
version of the data reported in this paper. All model calculations except WHDG [26] use a hydrodynam- 
ical description of the medium, but different extrapolation assumptions from RHIC to LHC. A variety 
of energy loss formalisms is used. An increase of Raa due to a decrease of the relative energy loss with 
increasing pj is seen for all the models. 

The curves labeled WHDG, ASW, and Higher Twist (HT) are based on analytical radiative energy loss 
formulations that include interference effects. Of those curves, the multiple soft gluon approximation 
(ASW (271) and the opacity expansion (WHDG [26]) show a larger suppression than seen in the mea- 
surement, while one of the HT curves (Chen [28]) with lower density provides a good description. The 
other HT (Majumder |[29l ) curve shows a stronger rise with pj than measured. The elastic energy loss 
model by Renk (elastic) QUI does not rise steeply enough with pj and overshoots the data at low px- 
The YaJEM-D model IlTTI . which is based on medium-induced virtuality increases in a parton shower, 
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shows too strong a pj -dependence of Raa due to a formation time cut-off. 

A more systematic study of the energy loss formalisms, preferably with the same model(s) for the 
medium density is needed to rule out or confirm the various effects. Deviations of the nuclear partem 
distribution functions (PDFs) from a simple scaling of the nucleon PDF with mass number A (e.g. shad- 
owing) are also expected to affect the nuclear modification factor. These effects are predicted to be small 
for pj > 10 GeV/c at the LHC [26] and will be quantified in future p-Pb measurements. 

4 Summary 

We have reported the measurements of charged particle pj spectra and nuclear modification factors 
Raa as a function of event centrality in Pb-Pb collisions at = 2.76 TeV. The results indicate a 

strong suppression of charged particle production in Pb-Pb collisions and a characteristic centrality and 
pj dependence of the nuclear modification factors. In central collisions (0 — 5%) the yield is strongly 
suppressed (Raa ~ 0.13) at pj = 6-7 GeV/c. Above pj = 7 GeV/c, there is a significant rise in the 
nuclear modification factor, which reaches Raa ~ 0.4 for pj > 30 GeV/c. This result is in agreement 
with the CMS measurement within statistical and systematic uncertainties. The suppression is weaker in 
peripheral collisions with Raa ~ 0.7 in the 70-80% centrality interval and is almost independent of pj. 
The observed suppression of high-px particles in central Pb-Pb collisions provides evidence for strong 
parton energy loss and a large medium density at the LHC. We observe that the suppression of charged 
particles with 5 < pj < 7 GeV/c reaches similar values when results from RHIC are compared to results 
from LHC in terms of the dA^h/drj. The measured Raa in 0-5% central collisions is compared to model 
calculations. An increase of Raa due to a decrease of the relative energy loss with increasing pj is seen 
for all the models. The measurement presented here, together with measurements of particle correlations 
ll32l and measurements using jet reconstruction [33], will help in understanding the mechanism of jet 
quenching and the properties of the medium produced in heavy-ion collisions. 
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